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Abstract 

In this paper we present a systematic analysis on the average interaction be- 
tween the last protons and neutrons in atomic nuclei, which can be extracted 
from the double differences of nuclear binding energies. The empirical aver- 
age proton-neutron interaction V^„ thus derived from experimental data can 
be described in a very simple form as the interplay of the nuclear mean field 
and the pairing interaction. It is found that the smooth behavior as well as 
the local fluctuations of the in even-even nuclei with N ^ Z are domi- 
nated by the contribution from the proton-neutron monopole interactions. A 
strong additional contribution from the isoscalar monopole interaction and 
isovector proton-neutron pairing interaction is seen in the Vpn for even-even 
N = Z nuclei and for the adjacent odd- A nuclei with one neutron or proton 
being subtracted. 

Keywords: Double binding energy differences, proton-neutron interaction, 
monopole, isovector pairing 

The binding energies of atomic nuclei reflect the interactions of its two 
constitutes, protons and neutrons. Differences of binding energies give nu- 
clear separation energies which can be used to isolate specific correlations. 
The zigzag behavior of one-body separation energies has long been well 
known. It provides clues to the pairing correlation between like nucleons [H, 
0]. In past decades, the structure of nuclei has been understood to a large 
extent within a mean-field (single-particle potential in the Hartree-Fock or 
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particle-hole channel) plus pairing approach. Modern nuclear structure model 
calculations within this framework can reproduce the bindingenergies of nu- 
clei over the whole nuclear chart with a high precision js], 0,1H. 

The correlation between the proton and the neutron has been expected 
to play a key role in the development of collective correlation joi 0] and in 
the evolution of the shell structure js], 0]. The (phenomenological) average 
interaction between the last protons and the last neutrons in even-even nuclei 
can be extracted from the double difference of binding energies as 

Vpn{Z,N) = \[B{Z,N) + B{Z-2,N-2) 



1 
4 

-B{Z-2,N)-B{Z,N-2)], (1) 



where B{Z, N) is the (positive) binding energy of a nucleus with Z pro- 
tons and N neutrons. The factor 1/4 takes into account the fact that four 
additional pairs are formed by the last two protons and neutrons. 

Recently, the proton-neutron interaction has attracted renewed interest 



which may reveal additional nuclear structure effects ll|, ll2|, |l3|, [ij, Il5|, 16 1 
and shed light on the possible existence of novel pairing correlation modes 171 . 

|21| . Stoitsov et al. showed that the global properties of Vpn can be 



reproduced by Hartree-Fock-Bogoliubov (HFB) calculations with the Skyrme 



functional plus a density-dependent 8 pairing interaction [1^ . A detailed caL 



culation was also done in Ref. [22| where the effects of the deformation and 
collective fluctuation on Vpn were analyzed. It would be interesting to un- 
derstand the microscopic mechanism behind the success of these calculations 
and to explore the extent to which the empirical Vpn can be incorporated into 
nuclear models where only proton-proton and neutron-neutron pairing corre- 
lations are explicitly taken into account. In this work I make an attempt in 
this direction by separating the contributions from the pairing and the mean 
field upon Vpn in a simple way. But perhaps even more appealing is to explore 
the local fluctuations of Vpn around the average values which large-scale HFB 



calculations fail to explain [ij] . These fluctuations may carry further nuclear 
structure information and serve as a constraint in future developments of 
nuclear structure models. 

The Vpn extracted from experimental nuclear binding energies Q (taken 
as positive values) are plotted in Fig. [H It can be seen that Vpn evolve rather 
smoothly as a function of mass number A. In fact, this average behavior of 
Vpn also probes the symmetry energy term (i.e., the isospin-dependence of 
the binding energy) in the macroscopic mass formula [24]. The overall trend 
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Figure 1: (color online). Empirical proton- neutron interactions in even-even nuclei ex- 
tracted from experimental nuclear masses [23| as a function of the mass number A. Open 
symbols correspond to calculations with the DZ [23] and HFB-17 [3] mass models. 



of Vpn can be well approximated by a smooth relation of {a + asA~^/^)/A (see 
the solid line in Fig. [1] that fits experimental data) [l^ . 

As a comparison, in Fig. [1] we also plotted the Vpn calculated from the 
HFB-17 [sl and Dufio-Zuker (DZ) jis], 26| mass models. We take the calcu- 
lated mass table from Ref. |3| for simplicity instead of repeating the large- 
scale HFB calculations done in Refs. 0, Q Q- The HFB-17 model is 
comprised of the conventional Skyrme functional, a contact pairing force and 
a number of empirical corrections including the Wigner energy. It can repro- 
duce experimental data within a deviation of 581 keV jsj. There are several 
versions of the DZ model available. Only the simplified DZ model is used 
in the present work. It contains ten terms and can reproduce experimental 
nuclear masses within a deviation factor of around 550 keV 0. It is seen from 
Fig. [Uthat both calculations can reproduce nicely experimental data, mostly 
within a deviation of 20%. 

We note that Eq. ([T]) can be rewritten as 

Wpn{Z,N) = S2n{Z,N)-S2niZ-2,N) 

= S2piZ,N)-S2piZ,N-2), (2) 
where S denotes the separation energy. From this relation it is easily seen 



^Onc may argue that the DZ model has more than ten parameters since it also contains 
several phenomenal scaling factors and isospin-dependent terms. 
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Figure 2: (color online). Empirical proton-proton (squares) and neutron-neutron (circles) 
interactions in even-even nuclei extracted from experimental nuclear masses as a function 
of the mass number A 



23j. The solid symbols denote those in the N = Z nuclei. 



that Vpn also measures the extra binding gained by the neutron (proton) pair 
when two additional protons (neutrons) are added. 

The two-nucleon separation energies in even-even nuclei can be written 

as 



S2niZ, N) = 2Sn{Z, iV - 1) + A„(Z, iV), (3) 

where 

A„(Z, N) = B{Z, N) + B{Z, N-2)- 2B{Z, N-1). (4) 

To a large extent A(Z, N) measures the pairing interaction between the last 
two like nucleons 0, 27|. The proton and neutron A{Z,N) extracted from 
experimental binding energies are plotted in Fig. [2l It is seen that the 
A(Z, A^) values are mostly within 1 — 5 MeV. They roughly follow a A~^^^ 
scaling (i.e., the solid line in the figure) but the trend is not smooth as that 
in Vpn- 

To understand the influence of the large fluctuations in the pairing inter- 
action on Vpn, we rewrite Eq. ([1]) as 

Vp,^{Z,N) = ^[Sn{Z,N~l)-Sn{Z-2,N-l)] 

+ ^[A^{Z,N)-An{Z-2,N)]. (5) 
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The quantities 5Sn{Z, N-l) = Sn{Z, N-l)-Sn{Z-2, N-1) and (5„(Z, A^) = 
An{Z, N) — An{Z, N — 2) measure the isospin dependences of the one-body 
separation energy (the mean-field) and pairing interaction, respectively. One 
can easily see that 6Sn{Z, N — 1) (and 6Sp{Z — 1,N)) also reveals the average 
proton-neutron interaction between the last proton pair and odd neutron as 



Vpr,{Z,N-l) 



hs^{z,N-i] 



(6) 



= -[B{Z,N -1) + B{Z -2,N -2) 
- B{Z -2, N-l)- B{Z, - 2)] . 

Contributions from the two basic ingredients 6S and 6 on the empiri- 
cal proton-neutron interaction Vpn can be extracted from experimental nu- 
clear masses. The results are plotted in Fig. |3l It is seen that Vpn is domi- 
nated by the contribution from 6S. The 5„ and 6p values are comparatively 
small, mostly within \S\ < 100 keV. This indicates that the empirical proton- 
neutron interaction can to a large extent be understood as a mean- field effect. 
This mechanism is supported by our calculations with the DZ and HFB-17 



models. It is also consistent with the observation of Ref. 14 that HFB 



calculations on Vpn are insensitive to the different choices of pairing forces, 
if the local fluctuations in the pairing interactions are negligible, it should 



be 



VpniZ, N) ^ VpniZ, VpniZ - 1, iV), 



(7) 



where Z and N are even numbers. This is indeed the case, as can be seen 
from Fig. |3l 

In Figs. [T]& [3] only nuclei with N Z are plotted. The extracted proton- 
neutron interactions Vpn in N = Z nuclei are noticeably larger than those 
of the adjacent N Z nuclei, as seen in Fig. HJ indicating that there is an 
additional binding in these self-conjugate nuclei. This is often described as 

Its origin has been intensively inve stig ated 



the Wigner effect [24, l29|, |30|, iM 



in the past decade in terms of proton-neutron pairing correlation [2J| and 



spin-isospin symmetry [32|. Ref. [2^ found that the Wigner effect can not 



be explained in term of J = 1 isoscalar neutron-proton pair correlations. 

To analyze this feature we consider a system with protons and n^, 
neutrons in a single- ?' shell. We assume that the two-body interaction obeys 



a simple form [33 



V = a + bti-t2 + GPo, 



(8) 
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Figure 3: (color online). Empirical proton- neutron interactions in odd-even (Red) and 
even-odd (Green) nuclei (i.e., SS) extracted from experimental nuclear masses [Ij, [28| . 
The open symbols correspond to Sp and (5„. The solid line describes the average behavior 
of Vpn in even-even N ^ Z nuclei. 
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Figure 4: (color online). Experimental Vpn values of even-even N ^ Z nuclei (filled 
circles) and the adjacent odd-odd (squares) and odd- A nuclei (triangles). The filled and 
open triangles correspond to systems with one nucleon subtracted from and added to the 
even-even nuclei, respectively. The solid line labeled 1* describes the average behavior of 
Vpn in even-even N Z nuclei from Fig. [T] 2* and 3* denotes its twice and three time 
values. 
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where Pq denotes the monopole pairing interaction. G is the corresponding 
(negative) couphng strength. The first two terms, which do not depend on 
the angular momentum J, define the "averaged" monopole interaction. The 
isovector and isoscalar channels of the monopole interaction are given by 



m;T=l 



a + 6/4, 



(9) 



and 



Vm;T=o = a- 36/4. (10) 
In usual shell-model Hamiltonians the values of Vm-T=i are around zero while 



those of Vm-T=o are strongly attractive (see, e.g., Refs. [3J, l35|), indicating 
that h should have a positive sign [sGj]. The J = two-body matrix element 
is given as {j'^\V\f)j=Q^T=i = a + 6/4 + (2j + 1)G. The total energy of the 
system can be written analytically as 



1^ , ^, 6 

£/ = en H — n[n — 1 H — 
2 ^ '2 



r(r+i)- — 



(11) 



+G 



n — V 



{A] + %-n-v)- T{T + 1) + s{s + 1) 



where e denotes the single-particle energy. The total number of nucleon pairs 
is n{n — l)/2 with n 
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T is the total isospin of the system, v 
and s denote the seniority and the reduced isospin. 

For the ground state of an even-even nucleus we have T = \nT, — ny\/2, 
s = and f = 0. For even-even nuclei with n-,^ ^ riy, we have 

^Vm-T=l + 2(K n;r=0 - Vm;T=l) ^ b 

4 



pn 



a. 



(12) 



The minus sign in the first term takes into account that the binding energy 
and Vpn are defined as positive in the present work. On the other hand, in 
the case of n^r = rii, (i.e., N = Z) , we have 



4:Vm-T=l + 3(Vm;T=0 ~ Kn;T= 



pn 



G 
~2 



G 
2 ■ 



(13) 



The difference between the Vpn in Eqs 



(Ha) and (dl is -{Vn^■,T=o - 
values would result in a sud- 
den kink in the one-body separation energy S when approaching the N = Z 



K„;T=i)/4 - G/2 or 6/4 - G/2. The large 



pn 
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Figure 5: (color online). The T — 1 (dashed line) and T = (dotted line) monopole 
interactions between a neutron pair and a proton pair in a single- j shell. We have a 
total number of four interacting pairs among which only one has T — 1. The isovector 
proton-neutron pairing matrix element also contributes in the T — 1 channel. 

line (c.f., Eq. ([6])), which can not be reproduced by usual mean-field calcula- 



A schematic picture is plotted in Fig. |5]to understand further the coupling 
of protons and neutrons in N = Z nuclei. For two proton- neutron pairs in a 
single- j shell, we have three T = interaction pairs and one T = 1 interaction 
pair. The T = 1 pairing matrix element also contributes to the total binding 
energy. For a system with two neutrons and two protons in different shells 
we have two interaction pairs for both the T = and 1 channels [32]. 

For a/ = j,T = l/2 system with three particles in a single-j shell, we 
have V = 1 and s = 1/2. The Vpn for such a nucleus can be expressed in the 
same form as above. The empirical relation of Eq. ([7]) still holds for these 
self-conjugate nuclei. In reality we have 



where Z takes even values. 

The empirical interactions between the odd proton and odd neutron in 
odd-odd nuclei can be extracted from binding energies in a way similar to 
those of even-even and odd-A systems. The ground state of odd-odd N = Z 
nuclei may carry isospin quantum numbers T = or 1. For the lowest T = 




VpniZ, Z) ^ Vpn{Z, Z-l)^ Vpn{Z - 1, Z) 



(14) 
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state one may extract the proton-neutron interaction as 



Vpn{Z-l.Z-l) 

= B{Z -1,Z-1) + B{Z -2,Z -2) 
-B{Z -1,Z -2)- B{Z -2,Z -1) 

= I - a. (15) 

The results are also plotted in Fig. SI 

Fig. m indicates that the Vpn in even-even N = Z nuclei and the adjacent 
odd-A nuclei with one less nucleon are roughly twice as large as those in 
neighboring N Z nuclei, while the Vpn in odd-odd N = Z nuclei are three 
times as large as the average values in N ^ Z nuclei. This may be understood 
from Eqs. f lT2|) . f lT3|) and f|T5|) by assuming that a ~ 0. The figure suggests 
that in reality b should be positive. In medium mass and heavy nuclei, it 
should also be much larger than the pairing strength G. This is consistent 



with the results of empirical shell model calculations [3J, |35| . The Vpn in very 
light even-even N = Z nuclei are much larger than the average behavior . 
This may be due to the additional binding gained from the enhanced pairing 
energy which scales like A~^^^. In the spin-isospin SU(4) symmetry limit, 
the Vpn of N = Z nuclei are four times larger than those for N ^ Z [32 



There was also no difference between Vpn in even-even and odd-odd 



nuclei 32 



We made no attempt to fit the strengths of the monopole interactions 
for single-j systems. In practice, for 0^5/2 nuclei with A ~ 20, we have 



Vm-T=o 5.5 MeV |3| and G 0.8 MeV. For O/7/2 nuclei with A ~ 50, 

we have Kn;T=o ~ —1.6 MeV f35|. These values reproduce reasonably the 
differences between Vpn in N = Z and N Z nuclei. One may expect 
that the "residual" proton-neutron interaction beyond the isoscalar monopole 
interaction also influence the This residual interaction, which is more 
complicated since it breaks the seniority symmetry [21], is not taken into 
account in the present work. 

On the first glance one may say that the isovector proton-neutron pairing 
also has significant infiuence on the Vpn for even-even N = Z nuclei. But 
in real cases the pairing interaction (and other matrix elements) is strongly 
modified by the monopole interaction. To evaluate the residual effect of the 
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pairing interaction we rewrite Eq. fllip as 

2a -G 



E 



en + 



-n(n — 1) 



(16) 



+ 



b-2G 



r(r+i)- — 



+(j + l)G(n - t;) + G 



— -t; + s(s + l) 
4 



from which it can be seen that it is the term (j + \)Gv that may result in 
an odd-even staggering in nuclear binding energies. This suggested that the 
residual pairing term in macroscopic mass formulas may be written as 



EpOc2 



(17) 



where v = 1 for odd- A nuclei and v = 2 for the T = \N — Z\/2 ground state 
of odd-odd nuclei. There should be no additional gain in pairing energy when 
crossing the N = Z line. This is consistent with the results in Fig. |2] where 
no noticeable difference is seen between the pairing energies oi N = Z and 
other nuclei. 

The DZ mass model is constructed starting from a shell-model monopole 
Hamiltonian as 

Hm = Hm + Hs + Hd, (18) 

where Hm is the macroscopic part including the symmetry energy (propor- 
tional to T{T + 1)) and the pairing energy. The microscopic spherical term 
Hs and the deformed term Hd take into account the residual three-body and 
four-body correlations between valence nucleons in the open shell. The ex- 
pectation value of the Hamiltonian Hg is calculated by assuming the normal 
filling scheme of nucleons. The deformed Hamiltonian Hd takes into account 
the effect of the promotion of valence nucleons to the next shell [i^. The 
DZ model also contains a phenomenological T{T — 1/2) term (referred to as 
the "Wigner" term in Ref. 2^) and a T/A correction to the pairing energy. 



These terms have limited influence on and the binding energies. One 
can obtain a slightly better agreement with experiments by refitting the DZ 
model parameters to the up-to-date mass table [1^ without these two terms 
and with the pairing term as Eq. f|T7|) . The deviation from experiment is 
a = 0.537 MeV. In this simplified model the only term that explicitly de- 
pends on isospin T is the sjTiimetry energy which induces the kink between 



Vpn for N 



Z and N ^ Z nuclei as seen in Fig. m 
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Figure 6: (color online). Experimental and calculated proton- neutron interactions Vpn in 
even-even Pb and Po isotopes. Experimental data are taken from Refs. 23 . 



We calculated the contributions from different terms of the DZ Hamil- 
tonian on Vpn- It is thus found that Hm reproduce the bulk properties of 
the Vpn while the many-body term Hs may lead to fluctuations around their 
mean values. 

Recently the Vpn in nuclei around = 126 and Z = 82 shell closures have 



been intensively investigated [ll|, |28|, |38[ . Experimental and calculated Vpn 
values in Pb and Po isotopes are plotted in Fig. [61 The Vpn values in nuclei 
below and above Z = 82 follow similar evolution patterns to those of Pb and 
Po isotopes, respectively. As can be seen from the figure, the evolution of Vpn 



shows a clear bifurcation pattern [ll[. That is nuclei below and above Z = 82 
evolve in two distinct ways as a function of neutron number N. Energy 
density functional calculations with the SkP Skyrme force cannot describe 



this behavior 1^. A better agreement was obtained for Pb isotopes in Ref. 



22| by taking into account the influence of beyond-mean-field correlations. 



Calculations with the DZ and HFB-17 models are also plotted in Fig. [6] for 
comparison. The DZ model shows a better agreement with experiments than 
other calculations. 

Calculations with the DZ model show that the Vpn values in this mass 
region are dominated by contributions from the monopole terms in Hm and 
Hs, as can be seen from Fig. [71 

A sudden drop at = 92 was noted in the Vpn values of Er isotopes, which 
was explained in terms of transitions from spherical to deformed shapes in 
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Figure 7: (color online). DZ model calculations (circle) on V^„ values of even-even Pb and 
Po isotopes and the contributions from Hg (triangle) and Hm (square) of the monopole 
Hamiltonian. 



Ref. . This is supported by calculations with the HFB-17 and DZ mass 
models. The results of DZ model calculations are plotted in the left panels of 
Fig. [HI in comparison with experimental data. The evolutions of Vpn values 
in Yb and Hf isotopes show similar patterns (in the latter case the location 
of the drop moves to = 94). In the right panels of the figure we plotted 
the contributions of the spherical and deformed monopole terms Hg and Hd. 
It can be easily seen that the sudden drop is a result of the competing effect 
of spherical and deformed monopole terms Hg and Hd- 

In summary, a systematic analysis on the empirical proton-neutron inter- 
action Vpn is done from a simple perspective by describing it as the interplay 
between the mean field and the pairing interaction. The results are also 
compared with those derived from existing nuclear energy density functional 
and monopole shell-model calculations. It is found that the bulk properties 
of the Vpn are dominated by the contribution from mean field which can be 
estimated empirically from the one-nucleon separation energies. The pairing 
energy plays a relatively minor role. For the same reason the Vpn for odd-A 
nuclei are close to those of the neighboring even-even nuclei. 

The Vpn for N = Z nuclei and the odd-A nuclei with one nucleon sub- 
tracted from the even-even N = Z ones are much larger than others. We 
analyzed this feature within a simple seniority model by including the isovec- 
tor and isoscalar monopole interactions and the monopole pairing interaction. 
It is thus found that in these nuclei there is a strong additional contribution 



12 



I ' I ' I ' I ' I ' I ' I 




I . I . I . I ' I . I ' I 



90 92 94 96 98 100 

N 



90 92 94 96 98 100 

N 



Figure 8: (color online). Left: experimental and calculated proton-neutron interactions 
Vpn in even-even Er (upper) and Yb (lower) isotopes. Right: Contributions from and 
Hd on Vpn values. 



from the isoscalar monopole interaction and isovector proton-neutron pair- 
ing interaction. As a result, the Vpn for even-even and odd-odd N = Z are 
roughly two and three times stronger than those for N Z nuclei. In the 
DZ model the cusp is induced by the T dependent term in the symmetry 
energy which is proportional to T(T+ 1). 

We also analyzed the local fluctuations in Vpn around their mean values. 
In the Dufio-Zuker mass model, these fluctuations are understood as the 
interplay of three-body- and four-body-like monopole interactions between 
valence nucleons in the open shell. As examples, we applied the DZ model to 
investigate the Vpn evolution pattern in nuclei around = 126 and A^ = 92, 
where intensive attempts were done within the Hartree-Fock approach. 

The author thanks R. Liotta and R. Wyss for stimulating discussions and 
R. Liotta for his reading of the manuscript. This work has been supported 
by the Swedish Research Council (VR) under grant No. 621-2010-4723. 
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